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Abstract. This note reviews the Peano-Baker series and its use to solve the general linear 
system of ODEs. The account is elementary and self-contained, and is meant as a pedagogic 
introduction to this approach, which is well known but usually treated as a folklore result 
or as a purely formal tool. Here, a simple convergence result is given, and two examples 
illustrate that the series can be used explicitly as well. 



1. Introduction 
Consider the inhomogeneous linear initial value (or Cauchy) problem 

(1) x = Ax + b, with x(t Q ) = x , 

on some interval I C R that contains to, with time-dependent quantities x, A and b. Here, 
A(t) denotes a family of matrices and b{t) a vector-valued function, both of dimension d over 
R (or over C). For simplicity, we assume A and b to be continuous on I. When LA(i), A(s)] = 
for all t,s £ I (where we use the notation [A, B] = AB — BA), the unique solution of (pQ) is 

(2) x(t) = exp( ! A{t) dr) (x + I exp(- I A(a) da) b(r) dr 



which simplifies to x(t) = exp(J* t A(t) dr^j x for the homogeneous case; see Cor. 2.41]. 
Also, when A(t) = A is constant, one has exp( J"* A{t) dr) = e^~*o) A as usual; we refer to [2, 
Sees. 11 and 12] for a general exposition. 

The solution of ([T]) is still unique when [^4(i), A(s)\ ^ for t ^ s, but the solution is then 
given by the Peano-Baker series (PBS); see [HI Sec. 16.5], Sec. 7.5] or [U Sec. 1.3] for 
background. This approach goes back to Peano (1888; see [H]), and was further developed 
by Baker (1905; see [3]). Baker's paper also contains a more detailed account of the history 
of this approach. The PBS is well-known in principle, but not widely known, and appears 
mainly in engineering textbooks; compare [HEJH] for examples. Quite frequently, it is used 
in a purely formal manner (without convergence considerations), and it is also sometimes 
claimed [T7] that it is of little practical use for the solution. However, as one can learn from 
the biographical sketches of Delone in {6], one should neither discard analytical tools nor 
believe repeated claims without evidence. In fact, when new problems surface, many (old) 
tools turn out to be more useful than expected. 

It is thus the purpose of this short note to summarise the elementary properties of the 

Peano-Baker series and to give a simple and self-contained account. In particular, we do 
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not restrict ourselves to the abstract fundamental system (or matrix solution), but discuss 
its explicit series expansion with convergence. The need for it came up in the context of a 
problem in mathematical population genetics; see [18] for details. We add two little examples 
to demonstrate that it can be used explicitly as well. 

2. Homogeneous case: Volterra equation and convergence 

Let 4> A (t]t Q ) denote the fundamental system of the homogeneous part of (P), which is also 
the flow of the linear system. It satisfies the Volterra integral equation 

(3) <Z> A (i;i ) = 1 + f A(T)<P A (T;t )dT, with* 4 (t ;t ) = l. 
By means of a formal Picard iteration, this leads to 

(4) $ A (t;t ) = l + 
where 

T n {t) 

For convenience, we set I (t) = 1. By construction, we then have the recursion 

(5) Z n+1 {t) = f A{T)X n (r)dT 

Jt 

for n € No- Clearly, I n (t ) = S n0 l, in line with <P A (t ;t ) = 1. Eq. ([!]) is known as the 
Peano-Baker series (PBS) in control theory [8, p. 598], see also Sec. IV. 5], or as the 
series expansion in the context of the time-ordered exponential in physics [UJ. The PBS was 
recently also extended to the setting of time scales [5] . 

Let us consider the individual terms of the PBS more closely. 

Lemma 1. When A is continuous on I, the matrix functions I n are differentiable and satisfy 
t n+1 (t) = A(t)Z n (t), for all t £ I and n E N . 

Proof. By definition ([5]), each matrix function X n is (component- wise) continuous, and hence 
so is the product AX n for every n. Applying the fundamental theorem of calculus separately 
to each component of the matrix equation ([5]) proves the lemma. □ 

Let us show next (in modern terminology) that the PBS is nicely convergent in our finite- 
dimensional setting; compare [T2J Sec. 16.5] or [TJ Sec. 2.11] for a slighlty different account, 
|16j for further results in this direction, and [IP} Ch. 4.3] for background on matrix Taylor 
series. We begin by establishing a link to the usual exponential series. In the one-dimensional 
case, this is [H Cor. 1.3.1]. 



A(r)dT+ A( Tl ) / A{t 2 )6.t 2 At x + ... = l + ^Z^i), 



n=l 



Mn) I A (r 2 ) ■ ■ l t M T n) dr„ • • • dr 2 dr x . 
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Lemma 2. If A is continuous on the interval I, with [A(t),A(s)] = for all t,s & I, one has 

Ut) = ^([ t A(r)dry 
Jt 

for all t £ I and n 6 N . In particular, this applies when A is one- dimensional. 

Proof. The claim is obviously true for n = and n = 1. We can now employ induction via 
the recursion ([5]), assuming the validity for some n € N. This gives 

l n+l {t) = f Air^Xnir^dT, = ! A(r 1 )^(/ Tl ^(r)dr) n dr 1 , 



"0 "0 



where the integrand on the right hand side can be rewritten as 
1 / /" r i \n 1 d / f T i 

-^)(/ /MdT ) .__(/^ (T)dT ) . 

This step employs the general chain rule, where the assumed commutativity is used. Inserting 
this expression into the previous formula completes the induction step by an application of 
the fundamental theorem of calculus. □ 

Remark 1. When [A(t), A(s)] = on /, the PBS (|4|) reduces to the well-known formula 

00 i ft ft 
*A(*;*o) =E^(/ ^(r)dr) =exp( / A(r) dr) 

n=0 ' 

as a consequence of Lemma El see [HI Sec. 2.3] for a detailed exposition of this case. This 
reduction is also mentioned in most of the sources cited so far; compare [U Cor. 1.3.2]. 

Remark 2. A closer look at the proof of Lemma [2] shows that the condition L4(t), A(s)] = 
may be replaced by the slightly weaker assumption that, for all t € I, the matrix A{t) 
commutes with the integral A(r) dr. 

Also, it should be mentioned that Lemma [2] can alternatively be proved by direct calcula- 
tions based on permutations of the integration variables followed by a suitable rearrangement 
to cover the integration region [t , t] n . When combined with an induction argument, it suffices 
to consider the permutations (12 . . . n), (1)(23 . . . n), . . . , (1)(2) . . . (n— 1, n). This approach is 
slightly more general. As it is also less transparent, we skip further details. 

Let now ||.|| denote any norm on W 1 (or on C d ), and define the compatible matrix norm 
as usual by \\A\\ := sup^n^ This implies the relations \\Ax\\ < ||j4|| • ||x|| and, more 

importantly, ||^4-B|| < ||j4|| • 

Theorem 1. If \\A(t)\\ is locally integrable on the interval I, the series representation of 
<P A (t;t Q ) is compactly convergent on I in the chosen matrix norm. 
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Proof. Let J C / be compact, with to € J. We show that the sequence of partial sums is 
Cauchy on J. So, let m, n G N with n > m and consider 

n m n n 

\\j2^(t) -J2 x ^t)\\ = || Yl ^ E INI 

k=0 k=0 k=m+l k=m+l 

/ / / "'iW-AMIdT^.dTadT!. 

k=m+l *0 *0 *0 

Since H-Afo) ■ ■ • -A(r fe ) || < H-A^))! • • • j4.(r fe ) , where all ||^4(Tj)|| are non-negative real num- 
bers, Lemma [2] implies that the last sum is majorised by 

£ t h(Jj A ^y, 

k=m+l l o 

which is the corresponding Cauchy estimate for the Taylor series of the ordinary exponential 
function on R, evaluated at J t * ||A(r)|| dr, which exists for all t S J by assumption. Since this 
series converges compactly, our claim follows. □ 

3. Solution of the inhomogeneous problem 

It is now obvious that <P A (t;t Q ) solves the homogeneous initial value problem This 
follows from a term-wise application of Lemma [T] to the PBS, which is fully justified by 
Theorem [TJ The determinant of <P A , which is a Wronskian and thus satisfies Liouville's 
theorem, reads 

(6) det($ A (t;t )) = det(#4(V,t )) -exp( I tv(A(r)) dr) = exp( j tr(^(r))dr), 

Jt Q Jt 

which never vanishes; see [21 Prop. 11.4] or [14} Thm. 2.23] for details. This means that <& A 
has full rank and thus indeed constitutes a fundamental system of the homogeneous linear 
system. 

Since <P A (t;t Q ) is the unique solution of ©, the flow property implies the relation 

(7) *A(f;*)*A(a;*o) = *A(t;t ), 

which, due to $ A (t;t) = 1, also implies (p A (t;s)) 1 = <& A (s;t). With the usual 'variation 
of constants' calculation, compare [21 Thm. 11.13], one can now easily derive the following 
result. 

Theorem 2. Let I be an interval, with t £ /. Suppose A is a continuous matrix function 
on I, and b is also continuous on I. Then, the Cauchy problem ([I]) has the unique solution 

ft 



x(t) = <P A (t; t ) (x + I $ A (t ; r) b(r) dr 



with <P A given by the PBS (jlj). When [A(t),A(s)] = for all t,s £ I, the formula simplifies 
to the explicit expression ([2]) with ordinary exponentials. □ 
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4. Examples 

Let us demonstrate the explicit applicability of the PBS with two simple examples. Both 
can also be solved by other means, but are perhaps still instructive. 

Consider the matrix family A(t) = ( J * ) with t E I, t = 0, and fixed parameter a. The 
matrices commute for a = 1, but not otherwise. One finds 







so that the PBS gives 



. . /V , . . e l - e at - (1 - a)te a 

*x(*;0) = L ;* I Wlth /(*)- 



e at I w (1 - a) 2 

Note that /(i) simplifies to ^i 2 e* for a = 1, in line with the then simpler ODE system. For 
general to, one determine <P A (t;to) from a similar calculation. The PBS differs both from the 
matrix exponential and from the Mathias formula for upper triangular matrices; compare [10, 
Thm. 3.6]. 

As a second example, consider the ODE system 

GK »)(:); 

which leads to the Airy function via the 2nd order ODE y = aty; see [TJ Ch. 10.4] for details. 
With a := a 1 / 3 , the PBS leads to the formula 



with 



9 (at) ^if(at) 
a 2 g(at) f(at) 



^y3/c 3 |~ 1 



/(*) =E 3fc (|) fc 7^T and 5 (z)=^3 fc (I) 



,_ w ^ to ^(u+iy. 

from [1, Eq. 10.4.3]. Here, (±) fe = |(± + 1) • • • (| + jfe - 1) as usual; see [TJ Eq. 6.1.22]. Note 
that 



*a(*;o) ^> ^ 

in line with the trivially solvable system for a = 0. 
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